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PROPOSAL OF A RISK MODEL FOR VEHICULAR TRAFFIC: A 
BOLTZMANN-TYPE KINETIC APPROACH 

PAOLO FREGUGLIA AND ANDREA TOSIN 


Abstract. This paper deals with a Boltzmann-type kinetic model describing 
the interplay between vehicle dynamics and safety aspects in vehicular traffic. 
Sticking to the idea that the macroscopic characteristics of traffic flow, includ¬ 
ing the distribution of the driving risk along a road, are ultimately generated 
by one-to-one interactions among drivers, the model links the personal (i.e., 
individual) risk to the changes of speeds of single vehicles and implements a 
probabilistic description of such microscopic interactions in a Boltzmann-type 
collisional operator. By means of suitable statistical moments of the kinetic 
distribution function, it is finally possible to recover macroscopic relationships 
between the average risk and the road congestion, which show an interesting 
and reasonable correlation with the well-known free and congested phases of 
the flow of vehicles. 


1. Introduction 

Road safety is a major issue in modern societies, especially in view of the con¬ 
stantly increasing motorisation levels across several EU and non-EU countries. Al¬ 
though recent studies suggest that this fact is actually correlated with a general 
decreasing trend of fatality rates, see e.g., [1, 32], the problem of assessing quant¬ 
itatively the risk in vehicular traffic, and of envisaging suitable countermeasures, 
remains of paramount importance. 

So far, road safety has been studied mainly by means of statistical models aimed 
at fitting the probability distribution of the fatality rates over time [24] or at fore¬ 
casting road accidents using time series [3, 21]. Efforts have also been made to¬ 
wards the construction of safety indicators, which should allow one to classify the 
safety performances of different roads and to compare, on such a basis, different 
countries [15, 16]. However, there is in general no agreement on which procedure, 
among several possible ones, is the most suited to construct a reliable indicator 
and, as a matter of fact, the position of a given country in the ranking turns out 
to be very sensitive to the indicator used. Despite this, a synthetic analysis is ulti¬ 
mately necessary: a mere comparison of the crash data of different countries may 
be misleading, therefore a more abstract and comprehensive concept of risk has to 
be formulated [28]. 

A recent report on road safety in New Zealand [2] introduces the following defin¬ 
itions of two types of risk: 

Collective risk: is a measure of the total number of fatal and serious injury 
crashes per kilometre over a section of road, cf. [2, p. 13]; 

Personal risk: is a measure of the danger to each individual using the state 
highway being assessed, cf. [2, p. 14]. 

While substantially qualitative and empirical, these definitions raise nevertheless an 
important conceptual point, namely the fact that the risk is intrinsically multiscale. 
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Each driver [microscopic scale) bears a certain personal level of danger, namely of 
potential risk, which, combined with the levels of danger of all other drivers, forms 
an emergent risk for the indistinct ensemble of road users [macroscopic scale). 
Hence the large-scale tangible manifestations of the road risk originate from small- 
scale, often unobservable, causes. Such an argument is further supported by some 
psychological theories of risk perception, among which probably the most popular 
one in the context of vehicular traffic is the so-called risk homeostasis theory. Ac¬ 
cording to this theory, each driver possesses a certain target level of personal risk, 
which s/he feels comfortable with; then, at every time s/he compares their per¬ 
ceived risk with such a target level, adjusting their behaviour so as to reduce the 
gap between the two [31]. Actually, the risk homeostasis theory is not widely accep¬ 
ted, some studies rejecting it on the basis of experimental evidences, see e.g., [11]. 
The main criticism is, in essence, that the aforesaid risk regulatory mechanism of 
the drivers (acting similarly to the thermal homeostatic system in warm-blooded 
animals, whence the name of the theory) is too elementary compared to the much 
richer variety of possible responses, to such an extent that some paradoxical con¬ 
sequences are produced. For instance, the number of traffic accidents per unit 
time would tend to be constant independently of possible safety countermeasures, 
because so tends to be the personal risk per unit time. Whether one accepts or 
not this theory, there is a common agreement on the fact that the background of 
all observable manifestations of the road risk is the individual behaviour of the 
drivers. In this respect, conceiving a mathematical model able to explore the link 
between small and large scale effects acquires both a theoretical and a practical 
interest. In fact, if on the one hand data collection is a useful practice in order 
to grasp the essential trends of the considered phenomenon, on the other hand the 
interpretation of the data themselves, with possibly the goal of making simulations 
and predictions, cannot rely simply on empirical observation. 

The mathematical literature offers nowadays a large variety of traffic models at 
all observation and representation scales, from the microscopic and kinetic to the 
macroscopic one, see e.g., [25] and references therein for a critical survey. Neverthe¬ 
less, there is a substantial lack of models dedicated to the joint simulation of traffic 
flow and safety issues. In [23] the authors propose a model, which is investigated 
analytically in [18] and then further improved in [22], for the simulation of car ac¬ 
cidents. The model is a macroscopic one based on the coupling of two second order 
traffic models, which are instantaneously defined on two disjoint adjacent portions 
of the road and which feature different traffic pressure laws accounting for more 
and less careful drivers. Car collisions are understood as the intersection of the 
trajectories of two vehicles driven by either type of driver. In particular, analytical 
conditions are provided, under which a collision occurs depending on the initial 
space and speed distributions of the vehicles. 

In this paper, instead of modelling physical collisions among cars, we are more 
interested in recovering the point of view based on the concept of risk discussed at 
the beginning. Sticking to the idea that observable traffic trends are ultimately de¬ 
termined by the individual behaviour of drivers, we adopt a Boltzmann-type kinetic 
approach focusing on binary interactions among drivers, which are responsible for 
both speed changes (through instantaneous acceleration, braking, overtaking) and, 
consequently, also for changes in the individual levels of risk. In practice, as usual 
in the kinetic theory approach, we consider the time-evolution of the statistical 
distribution of the microscopic states of the vehicles, taking into account that such 
states include also the personal risk of the drivers. Then, by extracting suitable 
mean quantities at equilibrium from the kinetic distribution function, we obtain 
some information about the macroscopic traffic trends, including the average risk 
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and the probability of accident as functions of the road congestion. To some extent, 
these can be regarded as measures of a potential collective risk useful to both road 
users and traffic governance authorities. 

In more detail, the paper is organised as follows. In Section 2 we present the 
Boltzmann-type kinetic model and specialise it to the case of a quantised space of 
microscopic states, given that the vehicle speed and the personal risk can be con¬ 
veniently understood as discrete variables organised in levels. In Section 3 we detail 
the modelling of the microscopic interactions among the vehicles, regarding them 
as stochastic jump processes on the discrete state space. Indeed, the aforemen¬ 
tioned variety of human responses suggests that a probabilistic approach is more 
appropriate at this stage. In Section 4 we perform a computational analysis of the 
model, which leads us to define the risk diagram of traffic parallelly to the more cel¬ 
ebrated fundamental and speed diagrams. By means of such a diagram we display 
the average risk as a function of the vehicle density and we take inspiration for pro¬ 
posing the definition of a safety criterion which discriminates between safety and 
risk regimes depending on the local traffic congestion. Interestingly enough, such 
regimes turn out to be correlated with the well-known phase transition between free 
and congested flow regimes also reproduced by our model. In Section 5 we draw 
some conclusions and briefly sketch research perspectives regarding the application 
of ideas similar to those developed in this paper to other systems of interacting 
particles prone to safety issues, for instance crowds. Finally, in Appendix A we 
develop a basic well-posedness and asymptotic theory of our kinetic model in meas¬ 
ure spaces (Wasserstein spaces), so as to ground the contents of the paper on solid 
mathematical bases. 

2. Boltzmann-type kinetic model with stochastic interactions 

In this section we introduce a model based on Boltzmann-type kinetic equa¬ 
tions, in which short-range interactions among drivers are modelled as stochastic 
transitions of microscopic state. This allows us to introduce the randomness of the 
human behaviour in the microscopic dynamics ruling the individual response to 
local traffic and safety conditions. 

In particular, we consider the following (dimensionless) microscopic states of 
the drivers: the speed u S [0, 1] C K and the personal risk u £ [0, I] C K, with 
u = 0, u = I standing for the lowest and highest risk, respectively. The kinetic 
(statistical) representation of the system is given by the (one particle) distribution 
function over the microscopic states, say / = f{t, u, u), t > 0 being the time 
variable, such that f{t, v, u) dv du is the fraction of vehicles which at time t travel 
at a speed comprised between v and v + dv with a personal risk between u and 
u -b du. Alternatively, if the distribution function / is thought of as normalised 
with respect to the total number of vehicles, /(t, v, u) dv du can be understood as 
the probability that a representative vehicle of the system possesses a microscopic 
state in [u, v + dv] x [u, u + du\ at time t. 

Remark 2.1. The numerical values of u introduced above are purely conventional: 
they are used to mathematise the concept of personal risk but neither refer nor imply 
actual physical ranges. Hence they serve mostly interpretative than strictly quant¬ 
itative purposes: for instance, the definition and identification of the macroscopic 
risk and safety regimes of traffic, cf. Section 4.2. The quantitative information on 
these regimes will be linked to a more standard and well-defined physical quantity, 
such as the vehicle density. 

Remark 2.2. The above statistical representation does not include the space co¬ 
ordinate among the variables which define the microscopic state of the vehicles. 
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This is because we are considering a simplified setting, in which the distribution 
of the vehicles along the road is supposed to be homogeneous. While being a 
rough physical approximation, this assumption nonetheless allows us to focus on 
the interaction dynamics among vehicles, which feature mainly speed variations, 
here further linked to variations of the personal risk. Hence the two microscopic 
variables introduced above are actually the most relevant ones for constructing a 
minimal mathematical model which describes trafhc dynamics as a result of con¬ 
current mechanical and behavioural effects. 


2.1. Use of the distribution function for computing observable quantities. 

If the distribution function / is known, several statistics over the microscopic states 
of the vehicles can be computed. Such statistics provide macroscopic observable 
quantities related to both the traffic conditions and the risk along the road. 

We recall in particular some of them, which will be useful in the sequel: 

• The density of vehicles at time t, denoted p = p{t), which is defined as the 
zeroth order moment of / with respect to both v and u: 


( 1 ) 


( 2 ) 


( 3 ) 


( 4 ) 


P(0 ■= / / u) dv du. 

Jo Jo 


10 JO 

Throughout the paper we will assume 0 < p < 1, where p = 1 repres¬ 
ents the maximum (dimensionless) density of vehicles that can be locally 
accommodated on the road. 

The average flux of vehicles at time t, denoted q = q{t), which is defined 
as the first order moment of the speed distribution, the latter being the 
marginal of / with respect to u. Hence: 




/(t, V, u) du I dv. 


The mean speed of vehicles at time t, denoted V = V{t), which is defined 
from the usual relationship between p and g, i.e., q = pV, whence 

Qit) 


ViJ) = 


pit)' 


The statistical distribution of the risk, say ip = ip(t, u), namely the marginal 
of / with respect to v: 

ip(t, u) := / f{t, V, u)dv, 

Jo 

which is such that (pit, u) du is the number of vehicles which, at time t, 
bear a personal risk comprised between u and u + du (regardless of their 
speed). Using ip we obtain the average risk, denoted U{t), along the road 
at time t as: 


■ [ uip(t, u) du = —r^ [ u( [ fit, V, u) dv] du. 

Jo pit) Jo \Jo J 

Notice that (p(t, u) du = p{t), which explains the coefhcient ^ in (4). 


W) 


2.2. Evolution equation for the distribution function. A mathematical mo¬ 
del consists in an evolution equation for the distribution function /, derived con¬ 
sistently with the principles of the kinetic theory of vehicular trafhc. 

In our spatially homogeneous setting, the time variation of the number of vehicles 
with speed v and personal risk u is only due to short-range interactions, which 
cause acceleration and braking. Since the latter depend ultimately on people’s 
driving style, they cannot be modelled by appealing straightforwardly to standard 
mechanical principles. Hence, in the following, speed and risk transitions will be 
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regarded as correlated stochastic processes. This way both the subjectivity of the 
human behaviour and the interplay between mechanical and behavioural effects will 
be taken into account. 

In formulas we write 


dtf = Q^{f, f)-fQ-{f), 


where: 

• Q^ifj /) is a bilinear gain operator, which counts the average number of 
interactions per unit time originating new vehicles with post-interaction 
state {v, u); 

• Q~{f) is a linear loss operator, which counts the average number of inter¬ 
actions per unit time causing vehicles with pre-interaction state (v, u) to 
change either the speed or the personal risk. 

Let us introduce the following compact notations: w := {v, u), I = [0, 1] C K. 

Then the expression of the gain operator is as follows (see e.g., [29]): 

(5) 

/)(i> w) = ri{w^, w*)T’(w* ^ w I w*, p)f{t, w*)/(t, w*) dw, dw* 

where w», w* are the pre-interaction states of the two vehicles which interact and: 

• 77(w*, w*) > 0 is the interaction rate, i.e., the frequency of interaction 
between a vehicle with microscopic state w* and one with microscopic state 
w*; 

• 7^(w, —>■ w|w*, p) > 0 is the transition probability distribution. More 
precisely, 'P(w* —>• w|w*, p) dw is the probability that the vehicle with 
microscopic state w* switches to a microscopic state contained in the ele¬ 
mentary volume dw of centred in w because of an interaction with the 
vehicle with microscopic state w*. Conditioning by the density p indicates 
that, as we will see later (cf. Section 3), binary interactions are influenced 
by the local macroscopic state of the traffic. 

For fixed pre-interaction states the following property holds: 


( 6 ) 


//2 


7^(w* —)• w I w*, p) dw =1 V w*, w* e Vp e [0, Ij. 


Likewise, the expression of the loss operator is as follows (see again [29]): 

Q~U){t,^)= [ p(w, w*)/(t, w*)dw*. 

Jp 

On the whole, the loss term fQ~{f) can be derived from the gain term (5) by 
assuming that the first vehicle already holds the state w and counting on average 
all interactions which, in the unit time, can make it switch to whatever else state. 
One has: 

f{t, w)Q"(/)(t, w) = r 7 (w, w*)'P(w ^ w' | w*, p)f{t, w)/(t, w*) dw' dw*, 

then property (6) gives the above expression for Q~{f). 

Putting together the terms introduced so far, and assuming, for the sake of 
simplicity, that 77(w*, w*) = 1 for all w*, w* G P, we finally obtain the following 
integro-differential equation for /: 


(7) 


dtf = // 'P(w* ^ w|w*, p)/(t, w*)/(t, w*)dw*dw*-p/. 
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Notice that, with constant interaction rates, the loss term is directly proportional 
to the vehicle density. Moreover, owing to property (6), it results 


f){t, w) - pf{t, dw = 0, 


112 

therefore integrating (7) with respect to w gives that such a density is actually 
constant in time (conservation of mass). 

For the study of basic qualitative properties of (7) the reader may refer to Ap¬ 
pendix A, where we tackle the well-posedness of the Cauchy problem associated 
with (7) in the frame of measure-valued differential equations in Wasserstein spaces. 
Such a theoretical setting, which is more abstract than the one usually considered 
in the literature for similar equations, see e.g., [6, Appendix A] and [26], is here 
motivated by the specialisation of the model that we are going to discuss in the 
next section. 

2.3. Discrete microscopic states. For practical reasons, it may be convenient 
to think of the microscopic states v, u as quantised (i.e., distributed over a set 
of discrete, rather than continuous, values). This is particularly meaningful for 
the personal risk u, which is a non-mechanical quantity naturally meant in levels, 
but may be reasonable also for the speed v, see e.g., [10, 12], considering that the 
cruise speed of a vehicle tends to be mostly piecewise constant in time, with rapid 
transitions from one speed level to another. 

In the state space /^ = [0, 1] x [0, 1] C we consider therefore a lattice of 
microscopic states {w^-}, with = {vi, Uj) and, say, i = 1, ..., n, j = 1, ..., m. 
For instance, if the lattice is uniformly spaced we have 


Vi = 


i — 1 
n — 1 ’ 


Ui = 


j - 1 


TO — 1 ’ 

and Vi < Vi+i for alH = 1, ..., n — 1, 


with in particular vi = ui = 0, Vn = Um = 1, 

Uj < Uj+i for all j = 1, ..., TO — 1. 

Proceeding at first in a formal fashion, over such a lattice we postulate the 
following form of the kinetic distribution function: 


( 8 ) 


/(i> w) = ^^/zi(t)(5w.,(w), 

i=l j=l 

where 5^^^ (w) = 6^^ (v) (g) Su^ (u) is the two-dimensional Dirac delta function, while 
fij{t) is the fraction of vehicles which, at time t, travel at speed Vi with personal 
risk Uj (or, depending on the interpretation given to /, it is the probability that a 
representative vehicle of the system possesses the microscopic state w^- = (vi, Uj) 
at time t). In order to specialise (7) to the kinetic distribution function (8), we 
rewrite it in weak form by multiplying by a test function (j) G C{ff) and integrating 
over P: 

^ ^ ^(w)/(t, w) dw 


dt 


'12 


= jj (()(w)7^(w*w I w*, p) dw^/(t, w*)/(t, w*) dw* dw* 

-p <('(w)/(t, w) dw; 

Jl'2 

next we read /(<, w) dw as an integration measure, not necessarily regular with 
respect to Lebesgue, and from (8) we get: 

n m 

i=l 
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= E E 

*.,**=1 


f/2 




D p)rfw ) 


i=l i = l 

In view of the quantisation of the state space, the transition probability distribution 
must have a structure comparable to (8), i.e., it must be a discrete probability 
distribution over the post-interaction state w. Hence we postulate: 




I Wi-j*, p) = 


i=l i=i 


where (p) S [0, 1] is the probability that the vehicle with microscopic state 

Wi.j, jumps to the microscopic state because of an interaction with the vehicle 
with microscopic state w^. , given the local traffic congestion p. Plugging this into 

the equation above yields 

n m 

i=i j=i 


= EE E E K 




^3 


whence finally, owing to the arbitrariness of </>, we obtain 


(9) 


•^4 = E E 


Remark 2.3. This discrete-state kinetic-type equation has been studied in the liter¬ 
ature, see e.g., [8, 10]. In particular, it has been proved to admit smooth solutions 
t I—>■ which are unique and nonnegative for prescribed nonnegative initial 

OTV/-1 in «->/-I/-I •! 1-1/-\n ■r'vT*/~>ir’/-i'*’T m I-T-m mricici \ \ ^ f. .( 

i=l J^3 \ 


data fij{0) and, in addition, preserve the total mass 'l2T=i in time. 


The arguments above can be made rigorous by appealing to the theory for (7) 
developed in Appendix A. In particular, we can state the following result: 

Theorem 2.4. Let the transition probability distribution have the form 

n m 

T’(w* ^ w| w*, p) = w*, p)d^.^{w), 

i=i j=i 

where the mapping (w*, w*, p) w*, p) is Lipschitz continuous for all 

i, j, i.e., there exists a constant Lip(7^*'^) > 0 such that 

|T’*'’(w*2, w^, P 2 ) W 2 , P2)\ 

< Lip(T’*-’)(^|w *2 - w*i| -t |W 2 - Wil -t |P 2 - /Oil) 

for all w,i, w, 2 , wj, W 2 € P, p € [0, 1]. 

Let moreover 

71 m 

/o(w) = E E (’«') 

i=ii=i 
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be a prescribed kinetic distribution function at time t = 0 over the lattice of micro¬ 
scopic states {wij} C P, such that 

n m 

f% > 0 Vi, j, 51XI = p e [0,1], 

j=ii=i 

SetVllj^ '■= p). Then the corresponding unique solution 

to (7) is (8) with coefficients fij{t) given by (9) along with the initial conditions 
fij{Q) = fij. In addition, it depends continuously on the initial datum as stated by 
Theorem A.5. 

Proof. The given transition probability distribution satisfies Assumption A.3, in 
fact 

Wi (T’(w*i ^ • I w*, Pi), 75(w*2 • I W 2 , P 2 )) 

n m 

= sup V Vp(Wy ) (T’*'’(w* 2, W^, P2) - W^, Pi)) 

<pGCi,,i(/=)nLipi(/2)““ 

( n m \ 

XX (|w* 2 -w*i| + |w;-w*| + |p 2 -pi|)- 

1=1 j=l ) 

Furthermore, /o G A4^(P). Then, owing to Theorem A.4, we can assert that the 
Cauchy problem associated with (7) admits a unique mild solution^. 

The calculations preceding this theorem show that if the fijitfs satisfy (9) 
then (8) is indeed such a solution, considering also that it is nonnegative (cf. Re¬ 
mark 2.3) and matches the initial condition /g. □ 

Theorem 2.4 requires the mapping (w*, w*, p) w*, p) to be Lipschitz 

continuous in PxPx [0, 1] but the solution (8) depends ultimately only on the val¬ 
ues p), cf. (9). Therefore, when constructing specific models, we 

can confine ourselves to specifying the values i*j*(p); taking for granted that 
they can be variously extended to points (w*, w*) ^ in a Lipschitz 

continuous way. 


3. Modelling microscopic interactions 

From now on, we will systematically refer to the discrete-state setting ruled 
by (9). In order to describe the interactions among the vehicles, it is necessary to 
model the transition probabilities (p) associated with the jump processes 

over the lattice of discrete microscopic states. 

As a first step, we propose the following factorisation: 

= Pi'oblMi. -t Uj I Tj., Vi*, p) ■ Prob(ui. -)• Vi \ Vi*, p) 

=: iV%,^,,{p)-(V")U*{p), 

which implies that changes in the personal risk (first term at the right-hand side) 
depend on the current speeds of the interacting pairs while the driving style (i.e., 
the way in which the speed changes, second term at the right-hand side) is not 
directly influenced by the current personal risk. In a sense, we are interpreting the 
change of personal risk as a function of the driving conditions, however linked to 
the subjectivity of the drivers and hence described in probability. By subjectivity 
we mean the fact that different drivers may not respond in the same way to the 
same conditions. More advanced models may account for a joint influence of speed 

^For the definition of mild solution to (7) we refer the reader to (15) in Appendix A. 
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and risk levels on binary interactions, but for the purposes of the present paper the 
approximation above appears to be satisfactory. 

As a second step, we detail the transition probabilities (p), (p) 

just introduced. It is worth stressing that they will be mainly inspired by a proto¬ 
typical analysis of the driving styles. In particular, they will be parameterised by 
the vehicle density p S [0, 1] so as to feed back the global traffic conditions to the 
local interaction rules. Other external objective factors which may affect the flow of 
vehicles and the personal risk, such as e.g., weather or road conditions (number of 
lanes, number of directions of travel, type of wearing course), will be summarised 
by a parameter a € [0, 1], whose low, resp. high, values stand for poor, resp. good, 
conditions. 

Remark 3.1. The interpretation of a is conceptually analogous to that of u dis¬ 
cussed in Remark 2.1: its numerical values do not refer to actual physical (meas¬ 
ured) ranges but serve to convey, in mathematical terms, the influence of external 
conditions on binary interactions. 

3.1. Risk transitions. In modelling the risk transition probability 

(^')ij..i-(p) = Prob(uj. ^ Uj I p) 

we consider two cases, depending on whether the vehicle with state {vi^, Uj^) in¬ 
teracts with a faster or a slower leading vehicle with speed Vi». 

• If < Vi* we set 

~ CiP^j,max{l, j, —1} T (1 ~ 

the symbol 6 denoting here the Kronecher’s delta. In practice, we assume 
that the interaction with a faster leading vehicle can reduce the personal 
risk with probability ap, which raises in high traffic congestion and good 
environmental conditions. The rationale is that the headway from a faster 
leading vehicle increases, which reduces the risk of collision especially when 
vehicles are packed (high p) or when speeds are presumably high (good 
environmental conditions, i.e., high a). Alternatively, after the interaction 
the personal risk remains the same with the complementary probability. 

• If Vi^ > Vi* we set 

("P m}) 

i.e., we assume that the interaction with a slower leading vehicle can only 
increase the personal risk because the headway is reduced or overtaking is 
induced (see below). 

3.2. Speed transitions. In modelling the speed transition probability 

(^")li*(p) = Pi'ob(u*. ^ v^ I Vi*, p) 
we refer to [27], where the following three cases are considered: 

• If Vi^ < Vi* then 

(^")lz*(p) = a(l - + (1 - a(l - p))5i,^,, 

i.e., the vehicle with speed Vi^ emulates the leading one with speed Vi* by 
accelerating to the next speed with probability 0(1 — p). This probability 
increases if environmental conditions are good and traffic is not too much 
congested. Otherwise, the speed remains unchanged with complementary 
probability. 
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• If Vi^ > Vi* then 

= a(l - + (1 - a(l - p))Si,z*, 

i.e., the vehicle with speed Vi_^ maintains its speed with probability a(l — p). 
The rationale is that if environmental conditions are good enough or traffic 
is sufficiently uncongested then it can overtake the slower leading vehicle 
with speed vt*. Otherwise, it is forced to slow down to the speed Vi* and 
to queue up, which happens with the complementary probability. 

• If Vi^ = Vi* then 

)i^i*(.P) 1} d” n(I P)^i,in\n{i^+1, n] 

+ (I - a - (I - 2a)p)Si^i^. 

In this case there are three possible outcomes of the interaction: if environ¬ 
mental conditions are poor and traffic is congested the vehicle with speed 
Vi^ slows down with probability (1 — a)p', if, instead, environmental con¬ 
ditions are good and traffic is light then it accelerates to the next speed 
(because e.g., it overtakes the leading vehicle) with probability q:( 1 — p); 
finally, it can also remain with its current speed with a probability which 
complements the sum of the previous two. 


4. Case studies 

4.1. Fundamental diagrams of traffic. Model (9) can be used to investigate the 
long-term macroscopic dynamics resulting from the small-scale interactions among 
vehicles discussed in the previous section. Such dynamics are summarised by the 
well-known fundamental and speed diagrams of traffic, see e.g., [20], which express 
the average fiux and mean speed of the vehicles at equilibrium, respectively, as 
functions of the vehicle density along the road. This information, typically obtained 
from experimental measurements [7, 19], is here studied at a theoretical level in 
order to discuss qualitatively the impact of the driving style on the macroscopically 
observable traffic trends. 

In Appendix A we give sufficient conditions for the existence, uniqueness, and 
global attractiveness of equilibria /oo G of (7), cf. Theorems A. 7, A. 8. 

Here we claim, in particular, that if the transition probability distribution V has 
the special form discussed in Theorem 2.4 then is actually a discrete-state 
distribution function. 


Theorem 4.1. Fix p G [0, 1] and let 7^(w* —>■ w | w*, p) be like in Theorem 2.4- 
Assume moreover that Lip(7^) < ^ (cf. Assumption A. 3). Then the unique equilib¬ 
rium distribution function f^o G which is also globally attractive, has the 

form /oo(w) = YJi=i Y.'T=i coefficients f°° satisfying 


( 10 ) 


XOQ 

Jij 


= - E 

p ^ 


E 


' 1 * 3 *- 




ip)fAufi 


u=l = l 


i = 1, 

j = 1, 


n 

, m. 


Proof. We consider directly the case p > 0, for p = 0 implies uniquely /oo = 0. 

Since Lip(7^) < i, we know from Theorems A.7, A.8 that (7) admits a unique 
and globally attractive equilibrium distribution /oo G A4(J.(/^), which is found as 
the fixed point of the mapping / i-A- ^Q^(/, /). In particular, defining the subset 

n m n m 

f G M(((/^) : /(w) = E E (w), fij > 0, ^ ^ fij = P 

2=1 j — 1 i—1 j — 1 



D := 
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it is easy to see that if V has the form indicated in Theorem 2.4 then the operator 
maps D into itself. In fact, for f G D we get 


1 

P 


«■"(/> /) = EE 

i=ii=i 



Y. Y. K,- 


.i* = li, J*=l 





i5w,,(w). 


From here we also deduce formally (10). Therefore, in order to get the thesis, it is 
sufficient to prove that D is closed in In fact this will imply that {D, Wi) 

is a complete metric space, and Banach contraction principle will then locate the 
fixed point of in D. 

Let (/^) C D be a convergent sequence in M^{P) with respect to the Wi 
metric. It is then Cauchy, hence given e > 0 we find gN such that ii h, k > 
then Wi (/^, /^) < e. This condition means 


112 


(/^(w) - /'"(w)) dw 


E E (4 - 4) 


i=i j=i 


< e 


for every (p G Cbp{l‘^) H Lip^)/^). In particular, taking a function ip which vanishes 
at every but one, say w^, we discover |</ 5 (wjj)||/jj — /A| < e for all h, k > N^. 
Thus we deduce that is a Cauchy sequence in K, hence for all i, j there exists 

fij G K such that —?> fij {k —>■ oo). Clearly > 0 because the f^Ps are all non¬ 
negative by assumption; moreover, Yh=i fij = linife-i-oo Yfi=i fij = P- 

Therefore / := YTi=i E^i ^ D. 

We now claim that /^ —>■ / in the Wi metric: 


W^i {f, f) = sup [ p{w) (/(w) - f’^iw)) dw 

veC'6,i(/2)nLipi(72)7z2 

n TTi 

sup EE(/b - fij) 

VGC6,i(/2)nLipi(72)^^ 


<EEI/*^-/' 

i=i j=i 


^ I k—^oo 


j 0 . 


This implies that D is closed and the proof is completed. 


□ 


Under the assumptions of Theorem 4.1, (10) defines a mapping [0, 1] 9 p i-A- 
{/“(p)}, i.e., for every p G [0, 1] there exist unique coefficients /“ solving (10) 
such that {/“} is the equilibrium of system (9) with moreover SjLi f^ — P- 
Owing to the argument above, for each p it is possible to compute the corres¬ 
ponding average flux and mean speed at equilibrium by means of formulas (1), (2) 
with the kinetic distribution fao- This generates the two mappings 

n m /X 

(11) P^ q{p) phAU(p):=^, 

i=l j=l p 

which are the theoretical definitions of the fundamental and speed diagrams, re¬ 
spectively, of traffic. Furthermore, it is possible to estimate the dispersion of the mi¬ 
croscopic speeds at equilibrium by computing the standard deviation of the speed: 


cry(p) : = 


\ 


- E -'^^p)f E 

P i=l j=l 


that of the flux being pavip), which gives a measure of the homogeneity of the 
driving styles of the drivers. 
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Figure 1 shows the diagrams (11), with the corresponding standard deviations, 
for different values of the constant a parameterising the transition probabilities, cf. 
Sections 3.1, 3.2. Each pair (p, q(p)), (p, V{p)) has been computed by integrating 
numerically (9) up to a sufficiently large final time, such that the equilibrium was 
reached. 

For a = 1 (best environmental conditions) the diagrams are the same as those 
studied analytically in [13]. In particular, they show a clear separation between the 
so-called free and eongested phases of traffic: the former, taking place at low density 
(p < 0.5), is characterised by the fact that vehicles travel at the maximum speed 
with zero standard deviation; the latter, taking place instead at high density (p > 
0.5), is characterised by a certain dispersion of the microscopic speeds with respect 
to their mean value. In [13] the critical value p = 0.5 has been associated with a 
supercritical bifurcation of the equilibria, thereby providing a precise mathematical 
characterisation of the phase transition. 

For a < 1, analogously detailed analytical results are not yet available and, to 
our knowledge. Theorems A.7, A.8 in Appendix A are the first results giving at 
least sufficient conditions for the qualitative characterisation of equilibrium solu¬ 
tions to (7) in the general case. According to the graphs in Figure 1, the model 
predicts lower and lower critical values for the density threshold triggering the phase 
transition. In addition to that, coherently with the experimental observations, cf. 
e.g., [19], some scattering of the diagrams appears also for low density, along with 
a capacity drop visible in the average flux (i.e., the fact that the maximum flux 
in the congested phase is lower than the maximum one in the free phase), which 
separates the free and congested phases as described in [33] . 

Compared to typical experimental data, the most realistic diagrams seem to be 
those obtained for a = 0.8, which denotes suboptimal though not excessively poor 
environmental conditions. It is worth stressing that such a realism of the theoretical 
diagrams is not only relevant for supporting the derivation of macroscopic kinematic 
features of the flow of vehicles at equilibrium out of microscopic interaction rules 
far from equilibrium. It constitutes also a reliable basis for interpreting, in a similar 
multiscale perspective, the link with risk and safety issues, for which synthetic and 
informative empirical data similar to the fundamental and speed diagrams are, to 
the authors knowledge, not currently available for direct comparison. 

4.2. Risk diagrams of traffic. Starting from the statistical distribution of the 
risk given in (3), we propose the following definition for the probability of accident 
along the road: 

Definition 4.2 (Probability of accident). Let u G (0, 1) be a risk threshold above 
which the personal risk for a representative vehicle is considered too high. We 
define the instantaneous probability of accident P = P(t) (associated with u) as 
the normalised number of vehicles whose personal risk is, at time t, greater than or 
equal to u: 

1 /.i ^ 

P{t) := - / y^{t, u) du = - / f{t, V, u)dvdu. 

P Ju P Ju Jo 

Asymptotically, using the discrete-state equilibrium distribution function f^ = 
foo{p) found in Theorem 4.1, we obtain the mapping 

1 ” 

(12) p^p{p)-.= - ^ 

j : Uj >11 i—1 

which shows that the probability of accident is, in the long run, a function of the 
traffic density. We call (12) the accident probability diagram. 
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P 



P 


Figure 1. Average flux (left column) and mean speed (right 
column) vs. traffic density for different levels of the quality of 
the environment. Red dashed lines are the respective standard 
deviations, n = 6 uniformly spaced speed classes have been used. 


Definition 4.2 and, in particular, (12) depend on the threshold u, which needs to 
be estimated in order for the model to serve quantitative purposes. If, for a given 
road, the empirical probability of accident is known (for instance, from time series 
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on the frequency of accidents, see e.g., [3, 21, 24]) then it is possible to find u by 
solving an inverse problem which leads the theoretical probability (12) to match 
the experimental one. This way, the road under consideration can be assigned the 
risk threshold u. 

The question then arises how to use the information provided by the risk threshold 
u for the assessment of safety standards. In fact, the personal risk u, albeit a prim¬ 
itive variable of the model, is not a quantity which can be really measured for each 
vehicle: a macroscopic synthesis is necessary. Quoting from [2]: 

Personal risk is most of interest to the public, as it shows the risk 
to road users, as individuals. 

To this purpose, we need to further post-process the statistical information brought 
by the kinetic model. Taking inspiration from the fundamental and speed diagrams 
of traffic discussed in Section 4.1, a conceivable approach is to link the personal 
risk, conveniently understood in an average sense, to the macroscopic traffic density 
along the road. For this we define: 

Definition 4.3 (Risk diagram). The risk diagram of traffic is the mapping (cf. (4)) 



(13) 


/oo = fooip) being the equilibrium kinetic distribution function of Theorem 4.1. 
The related standard deviation is 



Using the tools provided by Definition 4.3, we can finally fix a safety criterion 
which discriminates between safety and risk regimes of traffic depending on the 
traffic loads: 

Definition 4.4 (Safety criterion). Let u G (0, 1) be the risk threshold fixed by 
Definition 4.2. The safety regime of traffic along a given road corresponds to the 
traffic loads p G [0, 1] such that 


U{p) + <Ju{p) < u. 


The complementary regime, i.e., the one for which U{p) -\- erjj^p) > u, is the risk 
regime. 

Remark 4.5. An alternative, less precautionary, criterion might identify the safety 
regime with the traffic loads such that U{p) < u and the risk regime with those 
such that U{p) > u. 

Figure 2 shows the risk diagram (13) and the probability of accident (12) at 
equilibrium obtained numerically for various environmental conditions as described 
in Section 4.1, using the heuristic risk threshold u = 0.7. 

For a = 1 the average risk and the probability of accident are deterministically 
zero for all values of the traffic density in the free phase {p < 0.5). This is a 
consequence of the fact that, as shown by the corresponding diagrams in Figure 1, at 
low density in optimal environmental conditions vehicles virtually do not interact, 
all of them travelling undisturbed at the maximum speed. In practice, vehicles 
behave as if the road were empty and consequently the model predicts maximal 
safety with no possibility of collisions. In the congested phase [p > 0.5), instead, 
the average risk and the probability of accident rise suddenly to a positive value, 
following the emergence of the scattering of the microscopic speeds, see again the 
corresponding panels in Figure 1. Then they decrease monotonically to zero when 
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Figure 2. Average risk (left column), with related standard devi¬ 
ation (red dashed lines), and probability of accident (right column) 
vs. traffic density for different levels of the quality of the environ¬ 
ment. n = 6 and m = 3 uniformly spaced speed and risk classes, 
respectively, have been used 


p approaches 1, for in a full traffic jam vehicles do not move. Notice that the 
maximum of the average risk and of the probability of accident is in correspondence 
of the critical density value p = 0.5. 
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P b P 


Figure 3. a. Comparison of the curves (11), (12), (13) in the case 
a = 0.8. b. Determination of safety and risk regimes of traffic for 
a = 0.8, with corresponding probabilities of accident. 

For a < 1, the main features of the diagrams p i-A U[p) and p i-A P{p) described 
in the ideal prototypical case above remain unchanged. In particular, a comparison 
with Figure 1 shows that the maximum of both diagrams is still reached in corres¬ 
pondence of the density value triggering the transition from free to congested traffic, 
see also Figure 3a. This is clearly in good agreement with the intuition, indeed it 
identihes the phase transition as the most risky situation for drivers. It is worth 
remarking that such a macroscopically observable fact has not been postulated in 
the construction of the model but has emerged as a result of more elementary mi¬ 
croscopic interaction rules. For a = 0.6, 0.8, the average risk and the probability 
of accident take, in the free phase of traffic, a realistic nonzero value, which first 
increases before the phase transition and then decreases to zero in the congested 
phase. Notice that the standard deviation of the risk is higher in the free than in the 
congested phase, which again meets the intuition considering that at low density 
the movement of single vehicles is less constrained by the global flow. For a = 0.5, 
environmental conditions are so poor that U{p), P{p) —t I when p —>■ O’*", namely 
in correspondence of the maximum mean speed (cf. the corresponding panels of 
Figure 1). 

According to [2]: 

Personal risk is typically higher in more difficult terrains, where 
traffic volumes and road standards are often lower. 

By looking at the graphs in the first column of Figure 2, we see that the results 
of the model match qualitatively well this experimental observation; for low p and 
decreasing a the average personal risk tends indeed to increase. 

Again, the most realistic (namely suboptimal, though not excessively poor) scen¬ 
ario appears to be the one described by a = 0.8. In this case, cf. Figure 3b, the 
safety criterion of Dehnition 4.4, i.e., 

U{p) + (Tu{p) < 0.7, 

individuates two safety regimes for traffic loads p G [0, pi) and p G (p 2 , I], with 
Pi « 0.275 and p 2 « 0.51, respectively. The first one corresponds to a probability 
of accident P{p) < 27%, the second one to P{p) < 38%. Not surprisingly, the 
maximum admissible probability of accident in free flow (p < pi) is lower than the 
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one in congested flow {p > P 2 ), meaning that the safety criterion of Definition 4.4 
turns out to be more restrictive in the first than in the second case. This can 
be understood thinking of the fact that in free flow speeds are higher and the 
movement of vehicles is less constrained by the global flow, which imposes tighter 
safety standards. 


5. Conclusions and perspectives 

In this paper we have proposed a Boltzmann-type kinetic model which describes 
the influence of the driving style on the personal driving risk in terms of microscopic 
binary interactions among the vehicles. In particular, speed transitions due to en¬ 
counters with other vehicles, and the related changes of personal risk, are described 
in probability, thereby accounting for the interpersonal variability of the human 
behaviour, hence ultimately for the subjective component of the risk. Moreover, 
they are parameterised by the local density of vehicles along the road and by the 
environmental conditions (for instance, type of road, weather conditions), so as to 
include in the mathematical description also the objective component of the risk. 

By studying the equilibrium solutions of the model, we have defined two macro¬ 
scopic quantities of interest for the global assessment of the risk conditions, namely 
the risk diagram and the accident probability diagram. The former gives the average 
risk along the road and the latter the probability of accident both as functions of 
the density of vehicles, namely of the level of traffic congestion. These diagrams 
compare well with the celebrated fundamental and speed diagrams of traffic, also 
obtainable from the equilibrium solutions of our kinetic model, in that they predict 
the maximum risk across the phase transition from free to congested flow, when sev¬ 
eral perturbative phenomena are known to occur in the macroscopic hydrodynamic 
behaviour of traffic (such as e.g., capacity drop [33], scattering of speed and flux 
and appearance of a third phase of “synchronised flow” [19]). Moreover, within 
the free and congested regimes they are in good agreement with the experimental 
findings of accident data collection campaigns: for instance, they predict that the 
personal risk rises in light traffic and poor environmental conditions, coherently 
with what is stated e.g., in [2]. 

By using the aforesaid diagrams we have proposed the definition of a safety 
criterion., which, upon assigning to a given road a risk threshold based on the 
knowledge of real data on accidents, individuates safety and risk regimes depending 
on the volume of traffic. Once again, it turns out that the risk regime consists of 
a range of vehicle densities encompassing the critical one at which phase transition 
occurs. This type of information is perhaps more directly useful to the public 
than to traffic controlling authorities, because it shows the average risk that a 
representative road user is subject to. Nevertheless, by identifying traffic loads 
which may pose safety threats, it also indicates which densities should be preferably 
avoided along the road and when risk reducing measures should be activated. 

This work should be considered as a very first attempt to formalise, by a math¬ 
ematical model, the risk dynamics in vehicular traffic from the point of view of 
simulation and prediction rather than simply of statistical description. Several 
improvements and developments are of course possible, which can take advantage 
of some existing literature about kinetic models of vehicular traffic: for instance, 
one may address the spatially inhomogeneous problem [10, 12] to track “risk waves” 
along the road; or the problem on networks [14] to study the propagation of the risk 
on a set of interconnected roads; or even the impact of different types of vehicles, 
which form a “traffic mixture” [27], on the risk distribution. On the other hand, 
the ideas presented in this paper may constitute the basis for modelling risk and 
safety aspects also of other systems of interacting agents particularly interested by 
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such issues. It is the case of e.g., human crowds, for which a quite wide, though 
relatively recent, literature already exists (see [9, Chapter 4] for a survey), that in 
some cases [4] uses a kinetic formalism close to the one which inspired the present 
work. 


Appendix A. Basic theory of the kinetic model in Wasserstein spaces 


Equation (7), complemented with a suitable initial condition, produces the fol¬ 
lowing Cauchy problem: 


9t/= Q+(/,/)-p/, t>0, we/2 

/(O, w) = /o(w), w G /^ 


with the compatibility condition fj 2 /o(w) dw = p. Recall that P = [0, l]^ c 
is the space of the microscopic states. The problem can be rewritten in mild form 
by multiplying both sides of the equation by and integrating in time: 

(15) 

fit, w) = e"^‘/o(w) 

-I- J JJ ^(w*w I w*, p)/(s, w*)/(s, w*) dw* dw* ds, 

where we have used that, in view of (6), p = /j 2 fit, w) dw is constant in t. 

In order to allow for measure-valued kinetic distribution functions, as it happens 
in the model discussed from Section 2.3 onwards, an appropriate space in which to 
study (15) is X := (^([O, T]; M'^iP)), where T > 0 is a final time and M^iP) is 
the space of positive measures on P having mass p > 0. An element / G A is then 
a continuous mapping 1 1 —)■ fit), where, for all t G [0, T], fit) is a positive measure 
with /j .2 fit, w) dw = p. 

A is a complete metric space with the distance sup Wi if it), git)), where 

*e[o,T] 


(16) Wi if it), git)) = sup [ (p(w)(p(t, w) - fit, w)) dw. 

t^GLip^ Jj2 

is the 1-Wasserstein distance between fit), git) G Ad^iP). In particular, 
LiPi(^^) = {(P e CiP) : Lip((p) < 1}, 

Lip((p) denoting the Lipschitz constant of p. 


Remark A.l. The definition (16) of Wi follows from the Kantorovich-Rubinstein 
duality formula, see e.g., [5, Chapter 7]. However, since in Ad^iP) all measures 
carry the same mass and, furthermore, the domain P is bounded with diam(/2) < 2, 
the supremum at the right-hand side is actually the same as that computed over 
the smaller set Cb^iiP) C Lip;^(/2), where 

CbpiP) = {pGCiP) : ||^||^< 1 }, 
see [30, Chapter 1]. Hence we also have: 

(17) ITi if it), git)) = sup [ (p(w)(p(t, w) - fit, w)) dw. 

veC'6,i(/2)nLipi(72)d72 

In the following we will use both (16) and (17) interchangeably. 

Remark A.2. Let A > 0 and let p G CiP) with Lip(ip) < A. Then 

[ <p(w)( 5 (t, w) -/(t, w))dw = A / ^^(p(t, w) -/(t, w))dw 
J/2 7/2 A 

< AlTi if it), git)), 
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considering that ^ G Lip;^(/^). 

We will occasionally use this property in the proofs of the forthcoming theorems. 
Notice that, if A < 1, this estimate is stricter than that obtained by using directly 
the fact that ip G Lip 2 (/^) (i.e., the one without A at the right-hand side). 


To establish the next results, we will always assume that the transition probab¬ 
ility distribution V satisfies the following Lipschitz continuity property: 

Assumption A.3. Let 7^(w* —>■ • | w*, p) G for all w», w* G I^, p G [0, 1], 

where is the space of probability measures on P. We assume that there exists 

Lip(7^) > 0, which may depend on p (although we do not write such a dependence 
explicitly), such that 

Wi (T’(w*i ^ • I wj, p), T’(w ,2 ^ • I W 2 , p)) < Lip(T’) (|w *2 - w*i| -h |w 2 - w*|) 
for all w*i, w, 2 , w*, W 2 G P and all p G [0, 1]. 


A.l. Existence and uniqueness of the solution. Taking advantage of the mild 
formulation (15) of the problem, we apply Banach fixed-point theorem in X to 
prove: 

Theorem A.4. Fix p G [0, 1] and let fo G There exists a unique f G 

(^([O, -l-oo); At((.(/^)) which solves (15). 

Proof. We assume p > 0 for, if p = 0, the unique solution to (15) is clearly / = 0 
and we are done. We fix T > 0 and we introduce the operator defined on X as 

w) := e-^Vo(w) 


-b 






p)f{s, w*)/(s, w*) dw* dw* ds. 


Then we restate (15) as / = =A^(/), meaning that solutions to (15) are fixed points 
of on A. Now we claim that: 

• y'i.X) C A. 

Let / G A. The non-negativity of fo (by assumption) and that of V (by 
construction) give immediately that jX{f){t) is a positive measure for all 
t G [0, T], Moreover, a simple calculation using property (6) shows that 
the mass of S^{f){t) is 

[ =A^(/)(^, w) dw = pe“^‘-b p^ [ ds = p. 

Ji^ Jo 

Therefore we conclude that 5^{f){f) G for all t G [0, T], 

To check the continuity of the mapping 1 1 —>■ S^{f ){t) we define 

J^{ip){s)-.= JJ ip(w)7^(w*-> w| w*, p)dw^/(s, w*)/(s, w*)dw*dw*, 

for (f G Cbp{P) n Lip;^(/^), then we take ti, t 2 G [0, T] with, say, <1 < t 2 
and we compute: 


Wr{y{f)ih),.y{f){t2))= sup 




-Pt2 _ g-ptl 


IP 


<b’(w)/o(w) d-w 


-b f .f[p]{s) ds — f J^[(p](s) ds 

Jo Jo 


< p\e-P*^ - I 1 

< 3p^|t2 - til, 


p j ds^ -b p^e y ds 




20 


PAOLO FREGUGLIA AND ANDREA TOSIN 


where we have used the Lipschitz continuity of the exponential function 
and the fact that \^{(p){s)\ < p^. Finally, this says that S^{f) G X. 

• If T > 0 is sufficiently small then =5^ is a contraction on X. 

Let f,g£X,ipG Lipj(/^), and dehne 

(18) w*) := / (/ 3 (w) 7 ^(w* w I w*, p) dw. 

Notice preliminarily that, owing to Assumption A.3, 

|' 0 (w* 2 , W 2 ) - '(/’(w*!, W*)| < Lip(7^) (|w ,2 - W*i| + |W 2 - W*|). 

Then; 

f </’(w)(^( 5 )(t, w)-^(/)(t, w))dw 

= J JJ ■!/’(w*, w*)( 5 *p* - /*/*) dw* dw* ds 

(/*, f* being shorthand for /(s, w*), /(s, w*) and analogously p*, g*) 


= J gP(^ t) J ^ ^ '0(w,, w*)5* dw*^ (/ - 

+ J (d* 


/*) dw* ds 
— /*) dw* ds. 


Notice that 


W 2 ) - ' 0 (w*, w*))p* dw* 


< pLip('P)|w2 - Wil 


and that the same holds also for the mapping w* ha Jp V'(w*, w*)/* dw*. 
Hence we continue the previous calculation by appealing to Remark A.2 
(at the right-hand side) and to the arbitrariness of (p to discover: 


IFi (^(/)(t), < 2pUpiV) (/(s), ^(s)) ds 

^0 

< 2Lip(T’) (l-e"^‘) sup Wi g{t)), 

tG[0.T] 


whence hnally 

sup Wi(y{f){t), y{g){t)) <2Up(r) { 1 - sup Wi{f{t), g{t)). 

tG[0, T] te[0,T] 


From this inequality we see that: 

— if Lip {V) > 5 then it suffices to take T < ^ log 2 L^‘^-p^-i obtain that 
^ is a contraction on A; 

— if Lip {V) < 5 then =5^ is a contraction on X for every T > 0. 

Owing to the properties above, Banach fixed-point theorem implies the exist¬ 
ence of a unique fixed point f € X oi y which solves (15). If Lip (7^) < ^ this 
solution is global in time, whereas if Lip {V) > | it is only local. However, a simple 
continuation argument, based on taking f{T) G Al((.(/^) as new initial condition 
for f = T and repeating the procedure above, shows that we can extend it uniquely 
on the interval [T, 2T]. Proceeding in this way, we do the same on all subsequent 
intervals of the form \kT, {k + l)r], k = 2, 3, ..., and we obtain also in this case a 
global-in-time solution. □ 
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A.2. Continuous dependence. By comparing two solutions to (15) carrying the 
same mass p we can establish: 

Theorem A.5. Fix p G [0, 1] and two initial data /oi, /02 G Let /i, /2 G 

+ 00 ); be the corresponding solution to (15). Then: 

(fi(t), f2(t)) < (/oi, / 02 ) Vt > 0. 

Proof. ^ Let p G Lipj(/^). Using (15) we compute: 



p{w){f2{t, W) - /i(t, 


w)) dw <e (/oi, / 02 ) 


+ '0(w*, w*) (/ 2*/2 - /i*/r) dw* dw*, 

where w*) is the function (18) defined in the proof of Theorem A. 4. By 

means of analogous calculations we discover 


^1 (/i(^)) h{t)) < e (/oi, / 02 ) 

+ 2pLip(U) [ (/i(s), / 2 (s)) ds, 

Jo 

whence the thesis follows by applying Gronwall’s inequality. 


□ 


Remark A.6. From Theorem A.5 we infer that if Lip(U) < ^ then 

lim Wi{fi{t), f 2 {t)) =0. 

V + OO 

That is, all solutions to (15) approach asymptotically one another. This fact pre¬ 
ludes to the result about the equilibria of (15) proved in Theorem A.8 below. 


A.3. Asymptotic analysis. In this section we study the asymptotic trends of (7), 
in particular we give sufficient conditions for the existence, uniqueness, and attract¬ 
iveness of equilibria. It is worth stressing that equilibria of the kinetic model are at 
the basis of the computation of fundamental and risk diagrams of traffic discussed 
in Section 4. 

Besides the methods presented here, we refer the reader to [17] and references 
therein for other ways to study the trend towards equilibrium of space homogeneous 
kinetic traffic models and for the identification of exact or approximated steady 
states. 


A.3.1. Existence and uniqueness of equilibria. Equilibria of (7) are time-independent 
distribution functions /oo = /oo(w) G such that 

Q^ifoo, /oo) - Pfoo = 0. 

If p = 0 then it is clear that the unique equilibrium is the trivial distribution 
function f^o = 0. Assuming instead p > 0, from the previous equation we see that 
equilibria satisfy 

(19) /oo = -Q+(/oo, /oo), 

P 

i.e., they are fixed points of the mapping / 1 —>■ ^Q^{f, /). The next theorem gives 
a sufficient condition for their existence and uniqueness, relying on the Banach 
contraction principle in A4(j_(/^). 

Theorem A.7. Let Lip(U) < For all p G [0, 1], (7) admits a unique equilibrium 
distribution function /oo G 


^Throughout the proof, we will adopt the shorthand notations /* = f{t, w*), f* = f{t, w*). 
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Proof. Throughout the proof we will assume p > 0. 

The operator ^(5+ maps into itself, in fact, given / S it is 

clear that /) is a positive measure and moreover 

[ /)(w)dw = - //' /(w*)/(w*)dw*dw* = p. 

Ji^P pJJi^ 

Moreover we claim that, under the assumptions of the theorem, it is a contraction 
on Indeed, let f,gG M^{P) and <p S Lip^(/^), then: 

- [ </5(w) (Q+(5, p)(w) - Q+(/, /)(w)) dw 
P Ji^ 

= ~Ii^ (A (5(w*) - /(w*)) dw* 

+ ^ (A w*)f{w*)dw*'^ (5(w*) - /(w*)) dw,, 

where if is the function (18). From the proof of Theorem A.4, we know that both 
mappings w* i—>■ Jj 2 if{'w^, w*)p(w,) dw, and w, i—>■ /^2 w*)/(w*) dw* are 

Lipschitz continuous with Lipschitz constant bounded by pLip(7^), hence from the 
previous expression we deduce 

- [ V5(w) (Q+( 5 , p)(w) - Q+(/, /)(w)) dw < 2 Lip(d5)IFi (/, p). 

P d/2 

Taking the supremum over (p at the left-hand side yields finally 

Wi Qq+(/, /), 5)) < 2Lip(d^)ITi (/, p), 

which, in view of the hypothesis Lip(7^) < d, implies that is a contraction. 

Banach fixed-point theorem gives then the thesis. □ 

A.3.2. Attractiveness of equilibria. Under the same assumption of Theorem A.7, 
the equilibrium distribution function f^o is globally attractive. This means that all 
solutions to (14) converge to /oo asymptotically in time. The precise statement of 
the result is as follows: 

Theorem A.8. Let Lip(7^) < d. Any solution f £ C([0, -too); of (14) 

converges to /oo in the lUi metric when t —>■ +oo. 

Proof. Set /o(w) = /(O, w) £ Since /oo is the solution to (14) when the 

initial datum is /oo itself. Theorem A.5 implies 

Wi (f(t), /oo) < (/o, /^) 

and the thesis follows. □ 
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